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INTRODUCTION 


The  objective  of  this  program  is  to  provide  a  comprehensive  set 

of  measurements  of  electron,  negative-ion,  and  positive-ion  mobiJ  .ties 

together  with  attachment  and  ionization  rate  coefficients  for  dry  air 

and  mixtures  of  dry  air  with  water  vapor.  To  cover  the  range  of 

electric  fields  encountered  during  the  period  of  air  conductivity 

following  an  electromagnetic  pulse  event,  these  measurements  are  made 

for  values  of  E/N  (the  ratio  of  electric  field  to  gas  density)  ranging 

from  ~  10  Vm  to  ~  3  x  10  1  Vm  .  This  range  of  E/N  corresponds  to 

4  7 

values  of  electric  field  at  atmospheric  pressure  from  ~  10  to  ~  10  V/m 
and  covers  values  of  electron  mean  energy  from  near  thermal  to  values 
well  above  the  onset  of  Impact  ionization. 

The  specific  tasks  addressed  in  the  present  phase  of  the  program 

Include  completion  of  measurements  in  dry  air  and  documentation  of  all 

the  data  determined  for  dry  air.*  In  this  present  report  ve  describe 

work  carried  out  on  the  remaining  tasks  in  this  phase,  the  alignment  and 

calibration  of  the  optical  absorption  cell  for  monitoring  water  vapor 

concentrations  in  the  drift  tube,  and  the  measurement  of  swarm 

parameters  in  a  mixture  of  2Z  in  dry  air.  For  the  latter  task, 

—22 

measurements  have  been  carried  out  over  the  same  range,  6  *  10  <  E/N 

—19  2 

<  3  »  10  Vm  ,  as  that  covered  for  the  measurements  in  dry  air. 


*The  results  of  this  work  have  been  Issued  in  the  Air  Force  Weapons 
Laboratory  Technical  Report  No.  AFWL-TR-83-55. 
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II.  EXPERIMENTAL  APPROACH 


1 .  GENERAL  METHODOLOGY 

The  method  adopted  for  determining  the  electron  and  Ion  swarm 
parameters  In  humid  air  is  the  time-resolved  measurement  of  electron  and 
ion  waveforms  using  a  pulsed  drift  tube.  In  this  method,  a  pulse  of 
electrons  is  injected  from  the  cathode  into  a  region  of  uniform  electric 
field  E  containing  a  gas  of  known  density  N.  In  general,  a  distribution 
of  electrons,  negative  ions,  and  positive  ions  is  created  in  the  drift 
space  as  a  result  of  collisions  of  electrons  with  gas  molecules  during 
the  drift  time  of  electrons  to  the  anode.  The  electron  transport  and 
swarm  parameters  are  then  determined  directly  from  measurements  of  the 
time-resolved  currents  of  electrons  and  negative  ions  arriving  at  the 
anode  and  of  positive  ions  arriving  at  the  cathode.  Resolution  of  the 
arrival  spectra  of  each  charged  species  leads  to  che  most  direct  inter¬ 
pretation  of  the  measurements  and  is  the  approach  used  in  the  present 
apparatus. 

A  complete  description  of  the  drift  tube  and  of  the  determina¬ 
tion  of  swarm  parameters  from  the  measured  current  waveforms  together 
•fith  their  associated  measurement  uncertainties  has  already  been 
documented  (Ref.  1)  in  conjunction  with  the  measurements  in  dry  air. 
Paragraphs  2  and  3  describe  the  additional  apparatus  and  procedures 
required  for  acquiring  quantitative  data  in  humid  air. 

2.  IN-S1TU  MONITORING  OF  WATER  VAPOR  CONCENTRATION 

The  preparation  of  gas  mixtures  having  accurately  known  composi¬ 
tion  is  difficult  when  a  minor  constituent  of  the  mixture  is  a  polar 
molecule  (such  as  water  vapor)  because  of  physical  adsorption  of  the 
polar  molecule  to  surfaces  within  the  container.  Equilibrium  is 
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eventually  attained  when  the  adsorption  is  balanced  by  desorption  and  is 
usually  strongly  dependent  on  temperature.  In  the  case  of  gas  mixtures 
containing  water  vapor  at  room  temperature,  the  time  to  attain 
equilibrium  has  been  found  to  be  several  hours  (Ref.  2). 

Reproducible  gas  mixtures  containing  water  vapor  have  been 
prepared  by  first  allowing  water  vapor  into  the  mixing  chamber  to  a 
known  density,  freezing  out  the  water  vapor  in  a  cold  finger,  adding  the 
other  component  gas  to  a  known  density,  and  finally  allowing  the  cold 
finger  to  warm  up  to  release- the  water  vapor  (Ref.  2).  However,  tnere 
remains  soma  uncertainty  concerning  the  final  density  of  water  vapor  in 
the  mixture  due  to  possible  variations  in  the  adsorption  of  water  vapor 
to  surfaces  in  the  presence  of  the  other  components  of  the  mixture. 
Reproducibility  of  measured  swarm  parameters  as  a  function  of  time 
following  the  preparation  of  the  mixture  does  not  ensure  that 
replacement  of  adsorbed  water  vapor  by  other  molecules  in  the  mixture 
has  not  occurred  during  mixture  preparation. 

To  remove  this  ambiguity  from  the  present  measurements,  a  method 

has  been  developed  for  continual  in-situ  monitoring  of  the  water  vapor 

density  in  the  drift  tube  prior  to,  during,  and  following  the 

preparation  of  humid  air  mixtures.  The  method  involves  the  measurement 

of  the  optical  absorption  produced  In  a  beam  of  hydrogen  Lyman-ct 

radiation  traversing  a  well-defined  path  in  the  drift  tube  volume,  the 

absorption  being  proportional  to  the  density  of  water  vapor  along  the 

path.  The  choice  of  the  wavelength  of  the  radiation  is  governed  by  the 

requirement  that  the  absorption  by  water  vapor  be  much  larger  than  that 

by  the  other  constituents  of  the  mixture.  At  the  wavelength  121.56  nm, 

corresponding  to  that  of  hydrogen  Lyman-a,  the  absorption  cross  section 

—  1 7  2 

of  water  vapor  (rtHjO)  Is  (Ref.  3)  1.44  x  10  cm  (independent  of 

•20  2 

temperature)  whereas  that  of  oxygen  0(02)  is  (Ref.  4)  1.23  x  10  4Ucm  , 
and  nitrogen  Is  transparent.  In  the  present  application,  the  absorption 
path  length  is  defined  by  the  boundaries  of  an  absorption  cell  appended 
to  the  drift  tube  envelope. 
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a.  Absorption  Cell 

A  scale  drawing  of  the  absorption  cell  is  given  in  Figure  1 
where  the  numbers  refer  to  the  major  parts  listed  in  Table  1.  The  cell 
comprises  two  main  subassemblies:  (1)  the  variable  length  absorption 
volume  bounded  by  the  concave  mirror  (item  7)  and  the  vacuum  window 
(item  31),  and  (2)  the  optical  component  housing  (item  20).  The  entire 
assembly  is  mounted  on  a  4-1/2  in  diameter  Conflat  flange  (item  4)  which 
is  sealed  with  a  metal  gasket  to  a  corresponding  mating  flange  (item  4a) 
welded  in  the  drift-rtube  envelope.  Figure  2  shows  the  position  of  the 
absorption  region  relative  to  the  electron  drift  region. 

The  length  of  the  absorption  region  is  variable  from  0.1  to  5  cm 
by  moving  the  vacuum  window  ixially  with  respect  to  the  noncave  mirror. 
This  is  achieved  via  the  vacuum  sealing  bellows  (item  5)  which  is  welded 
at  one  end  to  the  inner  circumference  of  the  Conflat  flange  (item  4)  and 
at  the  other  end  to  the  cuter  circumference  of  a  custom  vacuum  flange 
(ir.em  15).  The  bellows  length  is  controlled  by  a  rotating  cylinder 
(item  14)  threaded  internally,  where  it  controls  the  position  of  the 
threaded  mating  drive  ring  (item  8)  which  is  screwed  firmly  to  the 
moving  vacuum  flange  (item  15).  Rotation  of  the  drive  ring  (item  8)  is 
prevented  by  three  guide  rods  only  one  of  which  (item  9)  is  shown  in 
Figure  1.  Rotation  of  the  drive  cylinder  (item  14)  is  possible  via  the 
ballbearing  race  (item  11)  chosen  to  accept  both  radial  and  axial  thrust 
and  for  its  high  precision.  The  outer  ring  of  this  ball  race,  and  the 
three  guide  rods  are  mounted  rigidly  on  the  support  flange  (item  10) 
which  is  rigidly  bolted  to  the  vacuum  flange  (item  4). 

For  vacuum  processing  of  the  drift  tube  the  optical  component 
housing  is  demounted  from  the  vacuum  flange  (item  15)  together  with  the 
bellows  drive  cylinder,  ballbearing  race  and  guide  rods  (items  14,  11 
and  9).  To  prevent  the  vacuum  window  from  collapsing  under  atmospheric 
pressure,  a  bakeout  fixture  is  attached  prior  to  the  removal  of  the 
drive  cylinder.  This  fixture  consists  of  four  bolts  (not  shown  in 
Figure  l)  which  are  fully  threaded  into  four  tapped  holes  in  the 
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Scale  drawing  showing  the  mechanical  details  of  the  single-ended  optical  absorption  cell 
The  numbers  used  to  label  the  various  parts  correspond  to  the  part  numbers  listed  in  the 
first  column  of  Table  1. 


TABLE  1.  LIST  OF  MAJOR  PARTS  OF  OPTICAL  ABSORPTION  CELL  SHOWN  IN  FIGURE  1 


P/N 

No. 

Used 

Description 

Supplier,  P/N 

4 

1 

4-1/2  Conflat  Flange 

a 

954-5079 

4a 

1 

A-l/2  Conflat  Flange  (Tapped) 

a 

954-5079 

5 

1 

Welded  Stainless  Steel  Bellows 

b 

622-60-3 

6 

1 

Welded  Stainless  Steel  Bellows 

b 

602-55-4 

7 

1 

2"  Diam.  VUV  Mirror  75  mm  FL 

c 

4436/7994 

8 

1 

Threaded  Drive  Ring 

9 

3 

1/4"  Diam.  Guide  Rods 

10 

1 

Support  Flange  for  Drive  Mech. 

11 

1 

Ballbearing  Race 

d 

KA-040-XP3 

14 

1 

Rotating  Drive  Tube-Intemally  Threaded 

15 

1 

Vacuum  Flange 

16 

1 

Gas-tight  Flange 

17 

1 

Gas-tight  Flange 

19 

3 

1/4"  Diam.  Support  Rods 

20 

1 

Optical  Cross  Housing 

21 

1 

Beam  Splitter  &  Lens  Holder 

24 

1 

O2  Filter  and  Spacer  Tube 

29 

2 

LiF  Biconvex  Lens ,  25  mm  FL 

e 

30 

1 

MgF2  Beam  Splitter  38  mm  Diam. 

e 

31 

1 

Tube-mounted  MgF2  Window 

e 

32 

2 

Photomultiplier  Housing 

f 

RFI/B-214  FV 

33 

1 

Lyman-a  Source  &  Exciter 

g 

399 

34 

2 

Solar-Blind  VUV  Photomultiplier 

f 

G26E  314  LF 

35 

m 

Insulating  Spacer  Flange 

36 

B 

Wedge  Flange 

37 

B 

Source  Aperture 

38 

B 

Image  Aperture 

a)  Varian  Corporation,  Palo  Alto,  CA 

b)  Metal  Bellows  Corporation,  Sharon,  MA 

c)  Oriel  Corporation,  Stamford,  CT 

d)  Keene  Corporation,  Kaydon  Bearing  Division,  Muskegon,  MI 

e)  Harshaw  Chemical  Company,  Solon,  OH 

f)  EMI  Gencom,  Plainview,  NY 

g)  Quantatec  International  Inc.,  Chatsworth,  CA 
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Scale  drawing  showing  the  location  of  the  absorption  region  between  the  mirror  and  the 
window  at  the  top  of  the  drawing  relative  to  the  electron  drift  region. 


drive  ring  (item  8)  and  whose  ends  bottom  on  the  support  flange 
(item  10)  to  maintain  a  spacing  between  the  retracted  window  and  the 
concave  mirror.  During  vacuum  processing,  the  mirror,  window,  bellows, 
drive  ring,  and  support  flange  (items  7,  31,  5,  8  and  10)  are  baked  at  a 
temperature  of  250°C  along  with  the  remainder  of  the  drift  tube  and 
associated  vacuum  system. 

The  optical  system  used  is  a  dual-beam,  single-ended  arrangement. 
Radiation  from  a  Lyman-o  source  (item  33)  is  split  into  two  paths  by  the 
beam  splitter  (item  30).  The  light  transmitted  by  the  beam  splitter  is 
focussed  at  the  center  of  curvature  of  the  concave  mirror  while  that 
reflected  from  the  front  surface  of  the  beam  splitter  is  focussed  at  the 
cathode  of  one  of  the  photomultiplier  tubes  (item  34,  upper).  The  light 
reflected  from  the  front  surface  of  the  beam  splitter  monitors  the 
intensity  emitted  from  the  source  and  is  used  as  the  reference  beam. 

The  light  transmitted  through  the  beam  splitter  passes  through  the 
vacuum  window  and  is  refocussed  by  the  concave  mirror  at  its  center  of 
curvature  and  thence,  after  reflection  from  the  back  surface  of  the  beam 
splitter,  at  the  cathode  of  the  other  photomultiplier  tube  (item  34, 
lower).  This  signal  beam  passes  through  the  absorption  region  twice,  so 
that  the  absorption  path  length  is  equal  to  twice  the  separation  between 
the  concave  mirror  and  the  vacuum  window.  The  focussing  properties  of 
the  system  are  kept  fixed  by  maintaining  a  constant  separation  between 
the  concave  mirror  and  the  remaining  optical  elements.  This  is 
accomplished  by  the  support  rods  (item  19)  which  rigidly  attach  the 
optical  cross  housing  (item  20)  to  the  support  flange  (item  10).  To 
accommodate  axial  motion  of  the  vacuum  window  within  the  constraints  of 
the  fixed  focus  system,  a  second  metal  bellows  assembly  (item  6)  is 
interposed  between  the  threaded  drive  ring  (item  8)  and  the  optical 
cross  housing  (item  20). 

To  avoid  unwanted  absorption  at  the  probing  wavelength  121.56 
nra,  the  optical  component  housing  external  to  the  drift  tube  is  made 
vacuum  tight  using  0-ring  seals  wherever  appropriate  and  evacuated  to  a 
pressure  ~  10  ^  Torr  using  a  separate  vacuum  system.  For  operation  at 
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121.56  nm  all  transmissive  optical  elements  (windows,  lenses,  beam 
splitter)  are  either  of  LiF  or  MgF2*  The  front  surface  of  the  concave 
mirror  is  aluminized  and  overcoated  with  MgF2«  The  beam  splitter,  which 
consists  of  a  MgF2  flat,  38  mm  diameter,  1  mm  thick,  is  mounted  in  a 
precision  machined  holder  (item  21)  which  also  carries  one  of  the  LiF 
lenses.  The  mounting  angle  of  45  deg  is  fixed  in  the  machining  of  the 
holder,  but  the  axial  position  and  rotation  of  the  holder  about  the  axis 
are  both  adjustable  within  the  optical  cross  housing  (item  20).  This 
adjustment  permits  accurate  alignment  of  both  beams  derived  from  the 
beam  splitter.  The  optical  cross  housing  is  machined  from  a  solid 
aluminum  block.  The  flanges  which  mate  with  the  photomultiplier 
housings  are  screwed  to  the  block  with  countersink  screws  and  a  vacuum 
seal  is  provided  by  the  incorporation  of  an  0-ring  in  each  joint. 

A  gas-tight  spacer  tube  (item  24),  joined  to  the  optical  cross 
housing  with  an  0-ring  seal,  couples  the  Lyman  a  source  (item  33)  and 
houses  the  source  aperture  (item  37)  and  the  collimating  lens.  The  1  mm 
diameter  source  aperture  provides  desirable  point-source  optical 
geometry  and  allows  the  use  of  focal  plane  apertures  for  discriminating 
against  scattered  light.  The  main  sources  of  scattered  light  in  the 
signal  beam  channel  have  been  identified  as  reflected  light  from  the  vacuum 
window  (item  31)  and  from  the  entrance  window  of  the  photomultiplier  which 
monitors  the  reference  beam  intensity.*  Scattered  light  from  the  vacuum 
window  is  reduced  considerably  by  the  4  mm  diameter  aperture  (item  38) 
placed  in  the  focal  plane  of  the  collimating  lens  and  coincidently  of 
the  concave  mirror.  Contributions  to  the  signal  beam  intensity  from 
light  reflected  by  the  reference  beam  photomultiplier  window  is  reduced 

*Reflection  also  occurs  from  the  entrance  -indow  of  the  photomultiplier 
which  monitors  the  signal  beam  intensity,  and  reflected  light  from 
this  source  contributes  to  the  reference  beam  signal.  However,  in  the 
present  system,  this  contribution  is  less  than  1Z  of  the  true  reference 
beam  intensity  even  with  no  absorption  in  the  drift  tube,  i.e.,  with 
maximum  signal  beam  intensity.  Consequently,  no  attempt  has  been  made 
to  reduce  this  scattered  light  component,  for  example,  by  tilting  the 
signal  beam  photomultiplier. 


by  tilting  the  photomultiplier  window  sufficiently  that  the  reflected 
light  is  outside  the  acceptance  angle  of  the  signal  beam  photomultiplier. 
This  is  accomplished  by  the  wedge  spacer  flange  (item  36)  interposed 
between  the  optical  cross  housing  (item  20)  and  the  photomultiplier 
housing  (item  32,  upper).  Both  photomultiplier  housings  are  joined  to 
the  optical  cross  housing  with  0-ring  seals.  In  addition  an  0-ring  seal 
between  the  entrance  window  of  each  photomultiplier  and  the  internal 
surface  of  its  housing  flange  preserve,  the  vacuum  integrity  of  the 
optical  component  housing. 

The  electronic  circuitry  associated  with  the  absorption  cell 
system  utilizes  photon  counting  and  is  illustrated  schematically  in 
Figure  3.  The  output  from  each  photomultiplier  is  fed  through  an 
amplifier  (LeCroy  Model  W100B)  having  a  voltage  gain  of  10  and  a  band 
width  of  230  MHz  to  a  discriminator  (LeCroy  Model  LD  601C)  having  a 
pulse-pair  resolution  of  9  ns  followed  by  a  320  MHz  decade  prescaler 
(Fairchild  Model  95H90)  and  then  to  one  channel  of  a  dual-channel 
counter.  The  counter  is  arranged  so  that  both  channels  remain  active 
until  a  predetermined  number  of  counts  have  been  accumulated  in  the 
reference  channel  B  at  which  time  the  counter  is  stopped  and  the 
contents  of  both  channels  read.  The  use  of  a  wide  bandwidth  system 
enables  a  statistical  variation  in  the  transmitted  count  of  -  1  part  in 

4 

10  to  be  obtained  with  a  counting  time  of  ~  10  s.  To  prevent 
electrical  ground  loops  between  the  absorption  cell  readout  circuitry 
and  the  waveform  recording  circuitry  from  the  drift  tube  an  insulating 
spacer  flange  (item  35  in  Figure  1)  is  inserted  between  the  optical 
cross  housing  (item  20  in  Figure  1)  and  the  bellows  flange  (item  17  in 
Figure  1). 

The  light  source  (item  33  in  Figure  1)  is  an  rf-excited  hydrogen 
discharge  lamp  which  emits  a  multiline  spectrum  together  with  an  exten¬ 
sive  continuum.  Two  filters  are  used  to  isolate  the  Lyman-a  line  from 
the  remaining  spectrum.  First,  the  photomultiplier  detectors  used  are 
of  the  solar-blind  type  and  are  sensitive  only  over  the  range  from  110 
nra  (the  MgF2  entrance  window  cut-off)  to  200  nm.  These  provide 
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broadband  filters  which  eliminate  radiation  of  wavelengths  longer  than 
200  nm.  Within  the  sensitive  range,  a  narrow  band  filter  transmitting 
at  the  Lyman-a  wavelength  is  provided  by  allowing  the  light  from  the 
source  to  pass  through  an  appropriate  optical  path  length  in  oxygen. 

The  gas-tight  spacer  tube  (item  24)  is  used  for  the  oxygen  absorption 
filter  cell.  For  convenience,  dry  air  is  substituted  for  pure  oxygen  in 
the  filter  cell  using  an  arrangement  in  which  the  cell  is  continuously 
flushed  at  atmospheric  pressure.  The  excellent  discrimination  provided 
by  this  filter  combination  is  illustrated  in  Figures  4  and  5.  Figure  4 
shows  the  spectrum  emitted  from  the  source  detected  by  the  solar-blind 
photomultiplier,  i.e.,  after  passage  through. the  broad  band  filter. 

This  spectrum  was  recorded  using  a  vacuum  monochromator  with  a  spectral 
resolution  of  0.3  nm.  To  accommodate  the  more  Intense  Lyman-a  line,  the 
photomultiplier  current  sensitivity  in  the  immediate  wavelength  vicinity 
of  Lyman-a  has  been  reduced  by  a  factor  of  10.  Figure  5  shows  the 
spectrum  recorded  by  the  photomultiplier  under  identical  conditions 
except  that  the  light  from  the  source  is  allowed  to  pass  through  a  dry 
air  filter  interposed  between  the  source  and  the  entrance  slit  of  the 
monochromator  and  having  the  same  optical  density  as  the  filter  cell  in 
Figure  1  (item  24).  Comparing  the  spectrum  of  Figure  5  with  that  of 
Figure  4  it  is  seen  that  while  the  intensity  of  the  Lyman-a  line  is 
reduced  by  only  ~  40%  the  remaining  unwanted  spectrum  is  reduced  to  the 
dark  current  level  of  the  photomultiplier.  Thus,  in  this  wavelength 
range  the  oxygen  filter  provides  a  much  narrower  bandwidth  filter  than 
is  available,  for  example,  from  a  dielectric  filter. 

The  modular  design  of  the  absorption  cell  system  enables  the 
complete  optical  alignment  of  the  system  to  be  carried  out  before 
attachment  to  the  drift  tube.  The  initial  alignment  is  carried  out 
using  the  light  emitted  from  the  source  which  lies  in  the  visible  region 
of  the  spectrum.  During  this  stage  the  LiF  lenses  are  replaced  by  glass 
lenses  whose  focal  length  in  the  visible  matches  that  of  the  LiF  lenses 
at  121.56  nm.  The  relative  position  of  the  beam  splitter  is  adjusted  so 
that  both  derived  beams  are  optically  centered  and  the  axial  position  of 
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Intensity.  Arbitrary  Units 


Figure  4.  Total  spectrum  emitted  from  the  Lyman-a  source  in  the  wavelength 
range  from  110  to  200  nm.  The  sensitivity  in  the  region  of  the 
Lyman-a  peak  is  reduced  by  a  factor  of  ten  compared  with  the 
remainder  of  the  scan. 
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Intensity.  Arbitrary  Units 
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Figure  5 
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.  Spectrum  recorded  under  identical  conditions  to  Figure  4  with  the 
exception  that  a  dry  air  filter  is  interposed  between  the  source 
and  the  detector.  The  sensitivity  in  the  region  of  the  Lyman— a 
peak  is  again  reduced  by  a  factor  of  ten.  The  remainder  of  the 
scan  represents  the  photomultiplier  dark  current. 


Che  image  aperture  adjusted  to  be  in  the  focal  plane  of  the  collimating 
lens.  In  the  final  alignment  stage  (with  the  LiF  lenses  installed)  the 
intensity  of  the  signal  beam  at  121.56  nm  is  optimized  by  adjusting  the 
lateral  position  of  the  source  with  respect  to  the  optic  axis  of  the 
cell  defined  by  the  source  aperture. 

The  variable  absorption  path  length  allows  differential  measure¬ 
ments  to  be  made  thereby  eliminating  end  effects  and  errors  associated 
with  condensation  of  water  vapor  on  the  vacuum  window  and  concave 
mirror.  In  addition,  by  making  differential  measurements,  the  true 
component  (i.e.,  that  component  which  has  passed  through  the  absorption 
region)  of  the  total  signal  detected  by  the  signal  beam  photomultiplier 
may  be  separated  from  the  unwanted  component  due  to  light  scattered  into 
the  detector  from  surfaces  within  the  optical  component  housing  as 
discussed  in  more  detail  in  paragraph  2b. 

b.  Calibration  of  Absorption  Cell 

For  a  beam  of  radiation  of  a  particular  wavelength  traversing  an 
optical  path  l  through  an  absorbing  gas  mixture  of  total  density  N,  the 
transmission  T  is  given  by 

T  -  exp  (-  l  p1Ni)  (1) 

where  and  are  the  absorption  cross  section  and  density,  respec¬ 
tively,  of  the  ith  component  species  of  the  mixture  and  N  *  pi^.  In  the 
present  application  of  absorption  of  Lyman-a  radiation  by  humid  air,  the 
absorbing  species  are  1^0  and  O2  and  Equation  (1)  may  be  written  in  the 
form 

T  -  exp  {-NZ  lf’(l-f)o(02)+  fa  (H2/0)]}  (2) 

where  f  is  the  fractional  density  of  water  vapor  in  the  mixture  and  f ' 
is  the  fractional  density  of  oxygen  in  dry  air,  i.e.,  f  «  [^Oj/N  and  f' 
■  (02 ] / [Dry  Air] . 
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Despite  precautions  taken  to  reduce  the  contribution  of  the 
measured  signal  due  to  scattered  light,  it  still  remains  a  substantial 
fraction  (~  40%)  of  the  total  signal.  However,  the  true  signal  com¬ 
ponent  which  has  passed  through  the  absorption  region  may  ba  extracted 
from  the  total  signal  by  making  differential  measurements.  If  I(£)  is 
the  total  intensity  of  the  signal  beam  detected  for  an  absorption  path 
l ,  then 


I(£)  -  Ii  +  I2  exp  (*  A  N  l)  (3) 

where  1^  is  the  intensity  contribution  due  to  scattered  light,  I2  is  the 
intensity  of  the  light  incident  on  the  absorption  region  in  the  drift 
tube,  and  A  is  given  [from  Equation  (2)]  by 

A  -  f'  (1-f)  o(02)  +  fo(H20)  (4 ) 

Equation  (3)  may  be  rewritten  in  the  form 

I(£)  -  I(£+A£)  exp  (ANA l)  -  I]_  [exp  (ANA£)-1 ]  (5) 

Thus,  by  making  measurements  of  the  total  intensity  of  the  signal  beam 
for  a  sequence  of  absorption  path  lengths  differing  by  a  constant 
increment  A£,  a  plot  of  I(£)  versus  I(i+A£)  is  linear  of  slope  S  given 
by 


S  -  exp  {NA£[f'(l-f)c(02)  +  fo(H20)]}  (6) 

Therefore,  provided  the  absorption  cross  sections  are  known,  the 
fractional  concentration  f  i  determined  from 

f  -  £nS/{NA£[o(H20)-f ’o(02)}  -  f 'o«>2)/ (o(H20)-f ' o«>2) ]  (7) 
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Using  published  values  of  Che  cross  sections  (Refs.  3,  4)  and  a 
value  of  f'  ■  0.209,  corresponding  Co  the  value  pertaining  to  the  dry 
air  used  in  the  present  study  (Ref.  1),  it  is  seen  that  the  ratio  of  the 
first  to  second  terms  in  Equation  (2)  is  less  than  12  for  values  of  f  * 
0.0173.  Thus,  for  the  present  application  in  humid  air  for  f  ■  0.02, 
Equations  (2),  (6)  and  (7)  may  be  replaced  with 

T  *  exp  (-Nifo(H20)]  (8) 

S  *  exp  [NAifo(H.20)]  (9) 

f  -  inS/[NAio(H20)]  (10) 

Calibration  of  the  present  absorption  cell  is  carried  out  with 

the  cell  attached  to  the  drift  tube.  The  procedure  consists  of 

introducing  a  known  density  of  pure  water  vapor  into  the  drift  tube, 

determined  using  the  Baratron  gauge  appended  to  the  system  (Ref.  1),  and 

measuring  the  intensity  of  Lyraan-a  radiation  transmitted  through  the 

absorption  cell  as  a  function  of  the  length  of  the  absorption  path. 

Typical  measurements  are  shown  by  the  solid  points  in  Figure  6  where  the 

intensity  1(1)  is  plotted  as  a  function  of  I(!+A!)  for  a  constant  path 

increment  A l  ■  0.8  cm  and  a  density  cf  water  vapor  N(H20)  »  1.421  x 

10  m  .  The  slope  of  the  plot  determined  by  linear  regression 

analysis  (as  indicated  by  the  solid  line  in  Figure  6)  gives  the  value  of 

a(H20)  from  Equation  (8).  From  a  series  of  similar  measurements  taken 

at  different  densities  of  pure  water  vapor,  the  cross  section  o(H20)  has 

been  determined  to  be  (1.42  t  0.3)  x  10  cm  ,  in  very  good  agreement 
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with  the  value  1.44  x  10  cm  obtained  in  previous  measurements  (Ref. 
3).  Using  this  value  of  o(H20)  the  fractional  concentration  of  water 
vapor  in  the  humid  air  mixture  Is  then  determined  from  Equation  (10). 


I  (i  +  A£) ,  Arbitrary  Units 


Figure  6.  Plot  of  che  intensity  of  the  Lyman-a  signal  measured  for  an 
absorption  path  SI  versus  that  measured  for  a  path  l  +  M. 

The  solid  points  correspond  to  measurements  in  pure  water 
vapor  and  the  open  points  to  measurements  in  humid  air 
(2.00%  H~0). 


'3.  PREPARATION  OF  DRY  AIR-WATER  VAPOR  MIXTURES 


The  gas-handling  and  vacuum  systems  associated  with  the  drift 
tube  facility  are  shown  schematically  in  Figure  7.  The  dry  air  and 
water  vapor  supplies  are  connected  to  the  gas  manifold  through  separate 
valves.  To  prepare  a  humid  ai  mixture,  the  minor  component,  i.e., 
water  vapor,  is  first  admitted  to  the  preevacuated  drift  tube  to  a  given 
pressure  and  the  drift  tube  isolating  valve  closed.  The  water  vapor 
pressure,  as  measured  with  the  Baratron  gauge  and  monitored  with  the 
optical  absorption  system,  decreases  initially  with  time  due  to 
adsorption  on  the  internal  surfaces  of  the  drift  tube.  After  a  period 
~  48  h,  the  pressure  finally  equilibrates  at  a  value  pw,  during  wh-'ch  time 
a  quantity  of  ~  10  Torr-2.  of  water  vapor  i3  adsorbed  within  the  drift 
tube.  The  gas  manifold  is  then  reevacuated  to  a  pressure  ~  10  ^  Torr  and 
pressurized  with  dry  air.  Dry  air  is  bled  into  the  drift  tube  through  the 
isolating  valve  to  a  total  pressure  p^.  given  by  pw/f,  where  f  is  the 
desired  fractional  concentration  of  water  vapor  in  the  mixture. 

After  preparation  or  the  mixture,  the  water  vapor  concentration 

is  monitored  using  the  absorption  cell.  Some  typical  measurements  are 

shown  by  the  open  points  in  Figure  6..  These  denote  the  intensity  I(£) 

plotted  versus  I(£+A£)  for  a  constant  path  increment  A£  »  0.8  cm  and  a 
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humid  air  mixture  of  total  density  N  *  71.0  x  10  m  formed  by  adding 

22  —3 

dry  air  to  a  water  vapor  sample  of  density  N(H20)  ■  1.421  x  10  m  , 
i.e.,  f  ■  0.0200.  Comparison  of  the  data  (solid  points)  taken  before 
the  addition  of  dry  air  with  the  measurements  (open  points)  recorded 
after  mixture  preparation  shows  that,  to  within  the  experimental 
uncertainty,  the  density  of  water  vapor  has  remained  unchanged.  This 
behavior  has  been  found  for  all  mixtures  prepared  in  the  present  study. 


The  residual  pressure  in  the  vacuum  and  gas-hagdling  systems  following 
bakeout  of  the  entire  system  is  in  the  low  10  Torr  range.  However, 
following  exposure  of  the  system  to  water  vapor,  the  base  vacuum 
pressure  attainable  without  rebaking  the  system  is  -  10  Torr  due 
to  slow  desorption  of  water  "apor  from  internal  surfaces. 


Moreover,  subsequent  periodic  absorption  measurements  have  ensured  that 
the  mixture  composition  remains  unchanged  during  measurements  of  swarm 
parameters  in  a  given  gas  mixture  sample. 

The  uncertainty  in  the  determination  of  the  fractional  concen¬ 
tration  of  water  vapor  from  the  optical  absorption  measurements  is 
estimated  to  be  ±  2%.  The  absorption  system  in  its  present  form  may  be 
used  for  the  determination  of  f  over  the  entire  range,  from  0.005  to 
0.04,  of  interest  for  the  measurement  and  swarm  parameters  in  humid  air. 
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III.  RESULTS  FOR  HUMID  AIR  (2Z  H20) 

The  present  measurements  have  been  carried  out  in  humid  air 

mixtures  with  f  *  0.0200.  The  dry  air  component,  supplied  by  Matheson 

Co.  is  obtained  from  the  same  cylinder  used  for  the  measurements  in  dry 

air  and  an  analysis  of  the  composition  has  been  given  previously 

(Ref.  1).  The  water  sample  was  prepared,  using  a  separate  ultrahigh 

vacuum  system,  by  vacuum  distillation  of  distilled  water  into  a  prebaked 

container  equipped  with  an  isolation  valve.  To  minimize  the  transfer  of 

impurities,  the  distillation  was  carried  out  near  room  temperature  and 

only  half  the  original  volume  was  transferred  into  the  container.  As  a 

further  check  on  the  purity  of  the  water  vapor  sample,  measurements  of 

electron  drift  velocity  in  water  vapor  have  been  carried  out  over  the 

—21  2 

range  of  E/N  from  10  to  50  x  10  V*m  .  The  present  measurements  are 
in  excellent  agreement,  to  within  the  combined  experimental  uncertain¬ 
ties,  with  the  currently  accepted  values  (Refs.  5,  6). 

1.  ELECTRON  MOBILITY 

A  summary  of  the  values  of  wg  and  the  derived  reduced  mobility 

Ue  as  a  function  of  E/N  is  given  in  Figure  8  and  in  Table  2.  The  upper 

abscissa  scales  in  Figure  8  (along  with  the  subsequent  figures)  show  the 

corresponding  values  of  electric  field  at  atmospheric  density  for  three 

different  altitudes,  0  km,  15  km,  and  30  km.  The  data  cover  the  range 
—22  —19  2 

of  E/N  from  6  x  10  to  3  *  10  V»m  .  The  uncertainty  in  the  values 

of  wg  and  ue  is  estimated  to  be  ±  IX  for  values  of  E/N  in  the  range 

10-20  <  E/N  <  10~19  V *m2  and  ±  2X  for  E/N  <  10'20  V*m2  and  E/N  > 

-19  2 

10  V*m  .  No  corrections  have  been  applied  to  the  data  for  attachment 
or  ionization  effects;  based  on  available  knowledge  such  corrections  are 
estimated  to  be  less  than  the  assigned  uncertainty. 
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TABLE  2.  VALUES  OF  THE  ELECTRON  NOBILITY  IN  HUMID  AIR  (2.00%  H2O) 
COMPARED  WITH  PREVIOUS  MEASUREMENTS  (REF.  1)  IN  DRY  AIR 


.  *  .  ■ 

wmm 

f(H20) 

»  0 

f  (h2o)  - 

0.0200 

E/N 

UeN 

yeN/N0 

UeN 

UeN/N0 

10-21v.m2 

1024/V-m*s 

m2/V‘s 

1024/V-m-s 

m2/V-s 

■■■■ 

j#v; .  ..r~  ^-ij 

1.613 

8.25 

7.43 

■ 

1.882 

7.59 

IS£9 

7.47 

1 

2.150 

7.10 

7.49 

1 

2.419 

6.78 

0.252 

7.56- 

2.688 

6.58 

0.245 

7.58 

SU 

1.250 

3.  "60 

6.04 

0.225 

7.41- 

0.276 

4.03 

5.39 

0.201 

7.40 

0.275 

- 

1.750 

4.70 

5.09 

0.189 

7.49 

0.279 

5.38 

4.77 

0.177 

7.38 

0.275 

■m 

6.72 

4.16 

0.155 

6.98 

0.260 

* 'U 

3.000 

8.06 

3.73 

0.139 

5.97 

0.222 

3.500 

9.41 

3.40 

0.127 

5.49 

0.204 

-  • 

4.00 

10.75 

3.18 

0.118 

5.05 

0.188 

4.50 

12.10 

2.96 

0.110 

4.51 

0.168 

i  ■  •' 

5.00 

13.44 

2.82 

0.105 

4.11 

0.153 

6.00 

16.13 

•2.58 

0.0961 

3.42 

0.127 

'  * 

7.00 

18.82 

2.39 

0.0887 

3.01 

0.112 

• 

8.00 

21.50 

2.24 

0.0832 

2.70 

0.100 

9.00 

24.19 

2.12 

0.0789 

2.49 

0.0926 

*.**•/' 

26.88 

2.09 

0.0777 

2.36 

0.0878 

\%  •  * 

12.50 

33.60 

1.94 

0.0723 

2.10 

0.0780 

40.3 

1.80 

0.0670 

1.96 

0.0729 

>■ 

17.  SO 

47.0 

1.71 

0.0638 

1  84 

0  0685 

■ 

53.8 

1.64 

0.0610 

1.78 

0.0662 

67.2 

1.56 

0.0582 

1.62 

0.0601 

V 

30.00 

80.6 

1.49 

0.0556 

1.48 

0.0549 

35.00 

94.1 

1.43 

0.0530 

1.43 

0.0531 

y" 

40.0 

107.5 

1.40 

0.0519 

1.39 

0.0518 

>ssm 

45.0 

121.0 

1.36 

0.0504 

H 

50.0 

134.4 

1.33 

0.0496 

1.33 

0.0496 

►  .  • 

60.0 

161.3 

1.24 

0.0461 

>‘v* 

70.0 

188.2 

1.21 

0.0449 

1.21 

0.0451 

-  /-■/ 

80.0 

215.0 

1.20 

0.0446 

90.0 

241.9" 

1.21 

0.0451 

f  *  '  ■ 

268.8 

1.21 

0.0450 

1.19 

0.0443 

336.0 

1.17 

0.0434 

403 

1.12 

0.0417 

1.13 

0.0419 

470 

1.12 

0.0417 

L- 

538 

1.10 

0.0407 

1.09 

0.0407 

f. 

f  .  •  .  * 

672 

1.01 

0.0375 

1.02 

0.0378 

1'  •'  ' 

• 

300.0 

806 

0.98 

0.0366 

0.98 

0.0365 

['•'  N  -  2.688  x  1025nf3 

f  0 


l  Estimated  uncertainties  in  the  measured  quantities  are: 

i  •'  E/N:  +0.2%;  f:  +2% 

V  .  •  WgN,  WeN/No  (f  =  0):  +  17.  for  E/N  _<  10“19  Vm2;  +  2%  for  E/N  >  10-19  V*m2 

JeN,  ueN/N0  (f  =  0.0200):  +  1%  for  10-20  _<  E/N  <_  10~19  V-ra2;  +  2%  for 

I;-/-  E/N  <  10-2°  V-m2  and  E/N  >  10"19  V-m2 
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A  comparison  of  the  present  values  of  electron  mobility  with 

previous  data  in  humid  air  for  f  ■  0.02  (Ref.  7)  and  in  dry  air 

(Refs.  1,  8)  is  presented  in  Figure  9.  For  the  cake  of  clarity  the  data 

in  dry  air  are  denoted  by  a  smooth  curve  drawn  through  the  experimental 

points.  The  left  ordinate  scale  represents  the  value  of  reduced  mobility 

while  the  right  ordinate  corresponds  to  the  mobility  at  a  density  NQ 

■  2.688  x  10^  m”^.  Also  indicated  in  Figure  9  (by  the  arrows  on  the 

left  ordinate)  are  the  values  of  the  mobility  of  thermal  electrons  in 

dry  air,  1.59  m  /V *s  deduced  from  the  measurements  of  electron  collision 

2 

frequencies  in  oxygen  (Ref.  9)  and  nitrogen  (Ref.  10),  and  1.75  m  /V»s 
deduced  from  measurements  of  electron  mobility  in  dry  air-carbon  dioxide 
mixtures  (Ref.  11). 

The  data  obtained  in  the  present  study  are  the  first  comprehen¬ 
sive  set  of  measurements  in  humid  air  over  an  extended  range  of  E/N. 

For  values  of  E/N  <  2  x  10”^  V»m^,  the  electron  mobility  is  found  to  be 

o 

constant  with  a  mean  value  of  0.278  ±  0.006  m  /V »s-  and  corresponds  to 
the  mobility  of  thermal  electrons.  This  value  is  in  very  good  agreement 
with  that,  0.271  ±  0.016  m^/V*s,  determined  by  van  Lint  (Ref.  7).  For 
values  of  E/N  >3  x  10  V»m  ,  the  data  agree,  within  experimental 
error  (±  2%),  with  the  previous  measurements  in  dry  air.  In  the 
intermediate  range  of  E/N,  the  values  of  van  Lint  (Ref.  7)  are  up  to  20% 
lower  than  the  present  results. 

As  may  be  seen  from  Figure  9,  the  effect  of  water  vapor  on  the 
electron  mobility  in  air  is  most  pronounced  at  low  values  of  E/N  where 
the  large  momentum  transfer  cr^ss  section  in  water  vapor  dominates 
despite  the  low  concentration  of  water  vapor  in  the  mixture.  The  major 
contribution  to  the  momentum  transfer  cross  section  at  very  low  electron 
mean  energies  is  expected  to  be  from  rotational  excitation. 


2.  NEGATIVE  ION  MOBILITIES 


Figure  10  where  Che  values  are  compared  with  previous  data  for  dry  air 
(Ref.  1).  As  in  the  case  of  dry  air,  the  negative-ion  waveforms  show 
the  presence  of  only  one  dominant  ion.  However,  based  on  the  widely 
different  mobilities,  the  identity  of  the  ions  for  values  of  E/N  <  A  x 
10  ^  V»m^  is  probably  different  from  that  for  E/N  >  6  x  10  ^V*m^  and 
the  dashed  line  in  Table  3  indicates  this  distinction. 

Over  the  range  of  E/N  from  6  x  10  ^  to  4  x  10  ^  V»m^,  the 
values  of  the  mobility  are  constant,  having  a  mean  value  of  (1.86  ± 
0.02)  x  10'4  m^/V*s.  This  suggests  that  over  this  range  of  E/N,  the 
observed  ions  are  in  thermal  equilibrium  with  the  gas.  The  lower  value 
found  for  the  mobility  in  humid  air  compared  with  that  in  dry  air 
indicates  that  the  dominant  ion  in  humid  air  has  a  different  identity 
than  that  in  dry  air.  In  this  region  of  E/N,  previous  work  (Ref.  12) 
has  shown  that  the  primary  ion  is  Of  formed  in  the  three-body  reaction 

e  +  02  +  M  *  Of  +  M  (11) 

where  the  third  body  M  is  either  an  oxygen,  nitrogen,  or  water  vapor 
molecule.  Among  the  possible  reactions  of  the  02  ions  which  can  occur 
are 


°2~  +  °2  +  M  *  °A~  +  M 

(12) 

Of  +  C02  +  M  ♦  C0f  +  M 

(13) 

Of  +  CO 2  ♦  C0f  +  02 

(1A) 

Of  +  H20  +  M  ♦  Of  •  H20  +  M 

(15) 

•  H20  +  H20  +  M  ♦  Of  (H20)2  +  M 

(16) 

0,"  •  H,0  +  CO,  ♦  C0A"  +  H,0 

(17) 

TABLE  3.  VALUES  OF  THE  NEGATIVE-ION  MOBILITY  IN  HUMID  AIR  (2.00%  H2O) 
COMPARED  WITH  PREVIOUS  MEASUREMENTS  (REF.  1)  IN  DRY  AIR 


E/N 

EN0/N 

f(H20)  =  0 

f  (h2o) 

=  0.0200 

U_N 

U-N/Nq 

y_N 

U-N/N0 

in-21_  2 
10  V«m 

104V/o 

lQ21/V-m-s 

10“4m2/V.s 

1021/V-m-s 

-4  2 

10  m  /V-s 

0.600 

1.613 

2.28 

4.97 

1.85 

0.800 

2.150 

2.27 

4.96 

1.85 

1.000 

2.688 

6.12 

2.28 

5.01 

1.86 

1.500 

4.03 

6.09 

2.27 

4.99 

1.86 

2.000 

5.38 

6.12 

2.28 

4.97 

1.85 

3.000 

8.06 

6.09 

2.27 

5.00 

1.86 

4.00 

10.75 

6.11 

2.27 

4.96 

1.85 

5.00 

13.44 

6.09 

2.26 

4.98 

1.85 

6.00 

16.13 

6.11 

2.27 

4.97 

1.85 

7.00 

18.82 

6.12 

2.28 

4.96 

1.85 

8.00 

21.50 

6.14 

2.28 

4.98 

1.85 

9.00 

24.19 

6.13 

2.28 

5.00 

1.86 

10.00 

26.88 

6.13 

2.28 

4.99 

1.86 

12.50 

33.60 

6.10 

2.27 

5.00 

1.86 

15.00 

40.3 

6.10 

2.27 

5.01 

1.86 

17.50 

47.0 

6.10 

2.27 

5.03 

1.87 

20.00 

53.8 

6.13 

2.28 

5.00 

1.86 

25.00 

67.2 

6.10 

2.27 

5.00 

1.86 

30.00 

80.6 

6.13 

2.28 

5.03 

1.87 

35.00 

94.1 

6.10 

2.27 

4.97 

1.85 

40.00 

107.5 

6.10 

2.27 

5.00 

1.86 

60.0 

161.3 

8.14 

3.03 

6.64 

2.47 

70.0 

188.2 

8.47 

3.15 

6.91 

2.57 

80.0 

215.0 

8.66 

3.22 

7.10 

2.64 

90.0 

241.9 

8.84 

3.29 

7.20 

2.68 

112.5 

302.4 

8.79 

3.27 

7.39 

2.75 

125.0 

336.0 

8.55 

3.18 

7.47 

2.78 

137.5 

369.6 

8.23 

3.06 

7.50 

2.79 

150.0 

403 

8.00 

2.98 

7.53 

2.80 

175.0 

470 

7.66 

2.85 

7.50 

2.79 

200.0 

538 

7.39 

2.75 

225.0 

605 

7.12 

2.65 

N  *  2.688  x  102^m  ^ 
o 

Estimated  uncertainties  in  the  measured  quantities  are: 
E/N:  +  0.2%;  f :  +  2% 

u_N,  u_N/Nq  (f  =  0):  +1% 

u_N,  u_N/N  (f  =  0.0200):  +  1% 


•M  lzoi 


(18) 


9 


02"  (H20)2  +  H20  +  M  ♦  02-  (H20)3  +  M  (18) 

Studies  of  these  reactions  have  been  made  with  O2  as  the  third 

body  (i.e.,  M-02 ) .  From  these  studies,  the  first  six  reactions  have 

thermal  rate  coefficients  of  3.0  x  10  43  m^/s,  4.7  x  10  4*  m^/s,  4.3  x 

10-16  m3/s,  2.2  x  10-40  m6/s ,  6  x  10"40  m6/s.  and  5.8  x  10“16  m3/s, 

respectively  (Refs.  13-16),  while  the  last  reaction  has  an  equilibrium 

3  -1 

constant  of  8.2  x  10  atm  (Ref.  14).  Assuming  that  the  same  values 
apply  with  N2  as  the  third  body  (i.e.,  M  ■  N2),  then  for  the  densities 


of  N2  (>  2.3  x  1024nf3),  02  (>  6.3  x  10‘JnfJ),  C02  «  3.6  x  1015m~J), 

22  —2 

and  H2C  (>  6  x  10  m  )  pertaining  to  the  air  mixtures  used  in  the 

•  _ 

present  study,  the  dominant  reaction  path  of  the  02  ions  is  expected  to 
lead  to  the  formation  of  02  (H^)^  via  the  reactions  described  by 

Equations  (i.5)  and  (16).  Further  clustering  of  this  ion  with  H20  is 
possible  as  described  by  Equation  (18).  The  forward  rate  for  this 
reaction  has  not  been  determined.  However,  from  the  value  of  the 
equilibrium  constant  for  this  reaction  and  the  range  of  water  vapor 
densities  given,  an  equilibrium  ratio  of  02~(H20)2/02~  (H20)2  >  18  is  to 
be  expected.  Thus,  if  the  forward  rate  for  02  (H20)^  formation  is 

comparable  to  that  for  02~(H20)2  formation  the  ion  species  observed  in 
the  present  work  is  expected  to  be  02  (^0)^.  Measurements  of  mass 

identified  ions  formed  in  mixtures  of  02  and  H20  confirm  this  conclusion 
(Ref.  12).  It  is  further  noted  that  for  humid  air  containing  naturally 
occurring  abundancies  of  C02  (~  0.03%)  and  fractional  concentrations  of 
H20  greater  than  0.2%  the  final  ion  identity  is  expected  to  remain 
unchanged. 

At  high  values  of  E/N,  the  primary  negative  ions  are  expected 
(Ref.  12)  to  be  0  and  H  formed  by  the  dissociative  attachment 
reactions 


23  -3  > 


,19  -3> 


e  +  02  +  0  +0 


e  +  H20  +  H_  +  OH 


I 


Further  reactions  involving  the  0  ions,  H  ions,  and  their  progeny  are 


0  +O2+M+O3  +  M  (21) 

0"  +  C02  +  M  ♦  C03"  +  M  (22) 

O3"  +  H20  +  M  ♦  O3-  •  H20  +  M  (23) 

03“  +  C02  ♦  C03'  +  02  (24) 

0“  +  H20  +  M  +  0“  •  H20  +  M  (25) 

H"  +  02  ♦  e  +  H02  (26) 

h"  +  H20  ♦  OH"  +  H2  (27) 

OH"  +  C02  +  M  ♦  OH"  •  C02  +  M  (28) 

OH"  +  H20  +  M  ♦  OH"  •  H20  +  M  (29) 


Studies  of  these  reactions  have  been  made  with  02  as  the  third  body 
(i.e.,  M  ■  02).  From  these  studies  the  thermal  rate  coefficients  have 
been  determined  to  be  1.1  x  10  43m^/s,  3.1  x  10  4^m^/s,  2.7  x  10  4®m^/s, 

5.5  x  10"l6m3/s,  1.3  x  10"40m6/s,  1.2  x  10"15m3/s,  3.7  x  10"l5m3/s, 

7.6  x  10  4®m^/s,  and  2.5  x  10  4®m^/s,  respectively  (Refs.  13,  14, 

17-19).  In  addition  to  the  above  product  ions,  more  complex  clusters  of 
0  ,  03”,  and  OH  with  H20  may  be  formed  (Refs.  12,  14). 

In  the  present  work  the  value  of  the  negative-ion  mobility  is 

—20  2 

found  to  be  a  function  of  E/N  for  values  of  E/N  >  6  x  10  V *m  .  This 
indicates  that  eitner  the  negative  ions  are  not  in  thermal  equilibrium 
with  tue  gas  and/or  that  the  identity  of  the  ion  species  is  a  function 
of  E/N.  It  is  also  noted  that,  based  on  the  lower  values  of  the 
mobility  in  humid  air  compared  with  those  in  dry  air,  the  presence  of 
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water  vapor  influences  the  ion  identity  in  this  range  of  E/N.  In  view 
of  the  liklihood  that  the  negative  ions  are  not  in  thermal  equilibrium 
with  the  gas  it  is  difficult  to  predict  with  certainty  the  identity  of 
the  ion  species  observed  in  the  present  work  based  on  arguments 
involving  thermal  reaction  rates.  Nevertheless,  if  these  rates  are 
taken  together  with  the  densities  of  Nj  (>  6  x  10  m  ),  (>  1.6  x 

10^m~^),  CO2  (<  6  x  101^m  J),  and  1^0  (>  1.6  x  lO^m  pertaining  to 
the  air  mixtures  used  in  this  range  of  E/N,  the  dominant  reaction  path 
of  ions  derived  from  0  is  expected  to  lead  to  the  formation  of  0  •  ^0 

and  higher  order  clusters  with  1^0  via  reactions  of  the  type  described 
by  Equation  ( 2 T ^ •  On  the  other  hand,  for  the  ratio  of  H2O/O2  ~  0.1 
pertaining  to  the  present  mixtures,  the  H-  ions  are  most  likely  to 
suffer  detachment  reactions  according  to  Equation  (26).  Thus,  even 
though  the  attachment  rate  coefficient  for  H-  formation  is  approximately 
a  factor  of  2  larger  than  that  for  0  formation  at  comparable  electron 
mean  energies,  the  contribution  of  H”  derived  ions  to  the  total 
negative-ion  spectrum  is  expected  to  be  small  in  the  present  mixtures. 

We  conclude  that  the  negative-ion  species  observed  in  the  present 
measurements  for  values  of  E/N  >  6  x  IQ  V*m  are  clusters  of  0  with 
H2O.  Cluster  formation  in  humid  air  has  a  significant  effect  on  the 
measurement  of  attachment  rates  and  is  discussed  further  in  paragraph  5. 

3.  POSITIVE  ION  MOBILITIES 

In  the  case  of  the  positive  ions,  the  waveforms  do  not  exhibit  a 
simple  exponential  time  dependence,  indicating  the  presence  of  more  than 
one  positive  ion  specie  or  the  occurrence  of  ion  conversion  reactions. 
Thus,  the  values  of  the  positive-ion  mobility  summarized  in  Figure  10 
and  Table  4  refer  to  the  slowest  ion,  the  only  specie  clearly  resolved 
in  the  present  measurements. 

The  positive  ions  formed  at  electron  energies  near  the  threshold 

X  x 

for  ion  production  by  electron  impact  are  expected  to  be  O2  ,  N?  ,  and 
as  a  result  of  the  reactions 
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TABLE  4.  VALUES  OF  THE  POSITIVE-ION  MOBILITY  IN  HUMID  AIR  (2.00%  H2O) 
COMPARED  WITH  PREVIOUS  MEASUREMENTS  (REF.  1)  IN  DRY  AIR 


E/N 

f (H20) 

*  0 

f (H20)  =»  0. 

O 

ro 

0 

0 

10“21V-m2 

104V/m 

d+n 

1021/V*m-  s 

y+N/N0 

-4  2 

10  m  /V-s 

y+N 

10237V*m-s 

y+N/No 

lO"4m2/V-s 

80.0 

215.0 

6.07 

2.26 

90.0 

241.9 

6.24 

2.32 

112.5 

302.4 

6.83 

2.54 

6.40 

2.38 

125.0 

336.0 

6.88 

2.56 

6.50 

2.42 

137.5 

369.6 

6.91 

2.57 

6.59 

2.45 

150.0 

403 

6.95 

2.59 

6.67 

2.48 

175.0 

470 

6.85 

2.55 

6.64 

2.47 

200.0 

538 

6.73 

2.50 

6.45 

2.40 

225.0 

605 

6.54 

2.43 

6.26 

2.33 

250.0 

672 

6.33 

2.35 

6.05 

2.25 

275.0 

739 

6.11 

2.27 

5.94 

2.21 

300.0 

806 

5.94 

2.21 

5.78 

2.15 

N  =>  2.688  x  1025m  3 
0 

Estimated  uncertainties  in  the  measured  quantities  are: 
E/N:  +  0.2%;  f  *  +  2% 
y+N,  U+N/Nq  (f  -  0):  +1% 

U ,N,  u . N/N  (f  =  0.0200):  +  1% 


(30) 


e  +  O2  *  02+  +  2e 
e  +  N2  ■*>  N2+  +  2e  (31) 

e  +  H20  +  H20+  +  2e  (32) 


The  more  probable  ion-molecule  reactions,  and  their  thermal  rate 


coefficients  (Refs.  16,  20-26),  that 
°2+  +  N2  +  N2  ♦  02+  •  N2  +  N2 

o2+  +  o2  +  o2  .-*•  o4+  +  o2 
o2+  •  n2  +  o2  *  o4+  +  o2 
o2+  +  ^0  +  n2  >  o2+  •  h2o  +  n2 
02+  +  H20  +  02  +  02+  •  H20  +  02 
°2+  •  N2  +  H20  -►  02+  •  H20  +  N2 
04+  +  H20  *  02+  •  H20  +  02 
02+  •  h2o  +  a2o  -►  h3o+  •  oh  +  o2 
n2+  +  n2  +  N2  ♦  N4+  +  N2 
n2+  +  02  >  02+  +  N2 

n2+  +  co2  >  co2+  +  n2 
co2+  +  o2  +  02+  +  co2 


can  occur  involving  these  ions  are 


(8  x 

10-43  m6/s) 

(33) 

(2.5 

x  10"42  m6/s) 

(34) 

(5  x 

10“17  m3/s) 

(35) 

(2.5 

x  10"40  m6/s) 

(36) 

(2.3 

c  10"40  m6/s) 

(37) 

(4  x 

10"15  m3/s) 

(38) 

(1.5 

x  10  m3/s) 

(39) 

(1.2 

x  10-15  m3/s ) 

(40) 

(5  x 

10"41)  m6/s) 

(41) 

(5.1 

x  10-17  m3/s) 

(42) 

(7.7 

x  10-16  m3/s) 

(43) 

(5.6 

x  10'17  m3/s) 

(44) 
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n4+  +  o2  ♦  02+  +  n2  +  N2 

(2.5 

x  10~16  m3/s) 

(45) 

n4+  +  fl20  +  h2o+  +  n2  +  n2 

(3  x 

10“15  m3/s) 

(46) 

N2+  +  H20  h20+  +  n2 

(2.2 

x  10  m3/s) 

(47) 

h2°+  +  o2  ■*>  o2+  +  h2o 

(2  x 

10-16  m3/s) 

(48) 

H20+  +  H20  +  H30+  +  OH 

(1.8 

x  10"15  m3/s) 

(49) 

H30+  +  H20  +  N2  ♦  h30+  •  H20  +  n2 

(3.4 

x  10“39  m6/s) 

(50) 

h3g+  +  H2O  +  o2  h3o+  •  h2o  +  o2 

(3.7 

x  10~39  m6/s) 

(51) 

H30+  •  OH  +  H20  ♦  H30+  •  H20  +  OH 

(1.4 

x  IQ'15  m3/s) 

(52) 

In  addition  to  these  product  ions,  more  complex  clusters  of  ^0*  with 
t^O  may  be  formed.  It  may  be  seen  that  these  reaction  schemes  for  the 
primary  ions  C^,  N2+,  and  H20+  all  lead  to  the  formation  of  clusters  of 
H30+  ions  with  water.  At  values  of  electron  energy  a  few  electron  volts 
above  the  thresholds  for  N2+,  02+,  and  H20+  formation,  N+,  0+,  and  0H+ 
ions  are  formed  by  dissociative  ionization.  Further  reactions  of  these 
ions  lead  to  the  additional  ion  species  N0+  and  C0+  and  clusters  of 
these  ions  with  H20. 

The  values  of  positive-ion  mobility  determined  in  humid  air  are 
found  to  be  a  function  of  E/N  and  also  are  different  from  the  values 
determined  in  dry  air.  Thus,  the  ion  identity  is  influenced  by  the 
presence  of  water  vapor  and  the  observed  ions  are  probably  not  in 
thermal  equilibrium  with  the  gas  over  the  range  of  E/N  covered  by  the 
present  measurements.  Thus,  while  the  application  of  thermal  reaction 
rates  to  the  interpretation  of  ion  identity  is  questionable,  the  common 
product  ion  of  the  different  reaction  schemes  makes  it  probable  that  the 
ions  observed  are  H30+  clusters  with  ri20. 
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4.  THREE-BODY  ATTACHMENT  COEFFICIENT 


Measurements  of  the  three-body  attachment  coefficient  have  been 

—21  2 

determined  over  the  range  0.6  <  E/N  <  30  x  10  V »m  and  a  summary  of 

these  values  is  given  in  Figure  11.  The  present  measurements  have  an 
estimated  uncertainty  of  ±  5%  and  the  values  reported  represent  the  mean 
of  three  or  more  measurements  at  each  value  of  E/N,  the  spread  being 
less  than  the  specified  uncertainty.  Measurements  carried  out  at 
different  gas  densities  for  a  given  value  of  E/N  have  established  the 
quadratic  dependence  of  the  attachment  coefficient  on  gas  density. 

In  Figure  11  the  present  measurements  are  compared  with  the  data 
of  van  Lint  (Ref.  7)  for  humid  air  having  the  same  fractional  concentra¬ 
tion  of  water  vapor  (2%)  and  with  the  previous  measurements  in  dry  air 
(Ref.  1).  The  values  of  the  effective  attachment  rate  coefficient 

summarized  in  Figure  12  and  in  Table  5,  are  obtained  by  combining  the 

2 

measurements  of  n/N-  with  the  measurements  of  wg.  These  values  have  an 
estimated  uncertainty  of  ±  IX. 

—21  2 

For  values  of  E/N  <2  x  10  V*b  the  attachment  rate 
coefficient  is  found  to  be  constant,  having  a  mean  value  of  (1.53  ± 

0.11)  x  10  m°/s.  Such  a  behavior  is  consistent  with  that  observed  in 

the  measurements  of  electron  mobility  and  further  corroborates  the 
interpretation  that  in  this  region  of  E/N,  the  electrons  are  in  thermal 
equilibrium  with  the  gas.  The  solid  arrow  on  the  left  ordinate  of 
Figure  12  indicates  the  value  of  the  thermal  attachment  frequency  in 
humid  air  (f  *  0.02)  deduced  using  recent  measurements  in  oxygen  and 
oxygen-nitrogen  mixtures  (Ref.  27)  and  in  oxygen-water  vapor  mixtures 
(Ref.  2),  i.e.,  from 

ka  "  Vn2  "  Wg/N2  -  0.0424  kj  +  0.159  k2  +  0.00412  k3  (53) 

Ir  Equation  (53),  kj,  k2  and  k3  are  the  thermal  attachment  rate 
coefficients  for  oxygen,  nitrogen  and  water  vapor,  respectively,  as  the 
third  body,  and  a  density  ratio  of  nitrogen  to  oxygen  to  water  vapor  of 
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Figure  11.  Comparison  of  the  present  measurements  of  the  three-body 
attachment  coefficient  as  a  function  of  E/N  with  previous 
data  in  humid  air  containing  the  same  water  vapor  concentration 
The  solid  curve  corresponds  to  previous  measurements  in 
dry  air. 


E  (NJN) ,  V/m  (N2  =  3.828x  10 


E  <  N./N) .  V/m  ( N,  =4. 049  x  10  m  .  216.7  K.  15  km) 


EIN./N).  V/m  IN.  =2.688  x  1Q~  m  273  K.  0  km) 


Humid  Air  (2%  H^O) 

•  Present 
7  van  Lint 


Figure  12. 


Comparison  of  the  present  measurements  of  the  three-body 
attachment  rate  coefficient  as  a  function  of  E/N  with  previous 
data  in  humid  air  containing  the  same  water  vapor  concentration 
The  solid  curve  corresponds  to  previous  measurements  in  dry  air 
The  dashed  and  solid  arrows  on  the  ordinate  correspond  to  the 
thermal  value  of  the  rate  coefficient  in  dry  air  and  humid  air 
(2%  H20) ,  respectively,  deduced  from  previous  measurements  in 
oxygen  and  mixtures  of  oxygen  with  nitrogen  and  water  vapor. 


TABLE  5.  VALUES  OF  THE  THREE-BODY  ATTACHMENT  RATE  IN  HUMID  AIR  (2.00% 
H20)  COMPARED  WITH  PREVIOUS  MEASUREMENTS  (REF.  1)  IN  DRY  AIR 


E/N 

10~21V*m2 

EN  /N 
o 

104V/m 

f (H20)  =  0 

f  (H20) 

-  0.0200 

1 

V  «>/■* 

SL 

. n-44  6, 

10  m  /s 

va(3)N02/N2 

107/s 

mam 

0.600 

1.613 

8.4 

11.7 

11.1 

15.3 

0.800 

2.150 

7.8 

10.7 

11.0 

15.3 

1.000 

2.688 

7.2 

9.9 

10.9 

15.1 

1.500 

4.03 

5.8 

8.1 

11.0 

15.2 

2.000 

5.38 

4.8 

6.6 

11.3 

15.7 

3.000 

8.06 

3.5 

4.8 

9.2 

12.8 

4.00 

10.75 

2.60 

3.6 

7.6 

10.5 

5.00 

13.44 

2.15 

2.97 

5.7 

7.9 

6.00 

16.13 

1.78 

2.46 

4.3 

5.9 

7.00 

18.82 

1.51 

2.09 

3.3 

4.5 

8.00 

21.50 

1.46 

2.02 

2.70 

3.7 

9.00 

24.19 

1.36 

1.88 

2.32 

3.2 

10.00 

26.88 

1.35 

1.87 

2.04 

2.83 

12.50 

33.60 

1.35 

1.87 

1.68 

2.33 

15.00 

40.3 

1.31 

1.81 

1.48 

2.05 

17.50 

47.0 

1.24 

1.71 

1.34 

1.86 

20.00 

53.8 

1.18 

1.63 

1.27 

1.75 

25.00 

67.2 

1.05 

1.45 

1.09 

1.51 

30.00 

80.6 

1.01 

1.40 

1.00 

1.39 

N  -  2.688  x  1025  nf3 
o 


Estimated  uncertainties  in  the  measured  quantities  are: 


E/N:  +  0.2%;  f:  +  2% 


'j  (3)/N2,  v.(3)N  2/N2(f  -  0):  +  67. 

a  a  O  “ 


.(3)/N2,  v  (3)N  2/N2  (f  =  0.0200):  +  7%. 
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77.4/20.6/2  has  been  assumed.  Thus,  using  the  values  of  k^  *  2.2  x 
10“42  m6/s,  k2  -  1.1  x  10"43  m6/s  (Ref.  27),  and  k3  -  1.4  x  10“41  m6/s 

"*  •  i  ^  /■ 

(Ref.  2),  gives  a  value  of  ka  ■  1.68  x  10  m°/s.  Also  shown  for 

comparison  in  Figure  12  by  the  dashed  arrow  on  the  left  ordinate  is  the 
thermal  attachment  frequency  in  dry  air,  1.15  x  10  m°/s,  deduced  from 

k  -  0.0441  k.  +  0.166  k0  (54) 

a  l  l 

using  the  same  values  of  k^  and  k2  and  assuming  a  density  ratio  of 
nitrogen  to  oxygen  of  79/21. 

The  present  value  for  the  thermal  attachment  frequency  in  humid 
air  is  30%  lower  than  that  determined  by  van  Lint  (Ref.  7)  but  is  in 
agreement,  to  within  the  combined  experimental  uncertainty,  with  the 
value  1.68  x  10”  m°/s  deduced  from  Equation  (53).  The  present  values 
of  attachment  frequency  at  higher  values  of.E/N  are  also  lower  than 
those  determined  by  van  Lint.  From  the  measurements  of  Hegerberg  and 
Crompton  (Ref.  27),  it  is  unlikely  that  the  attachment  cooling  effect 
may  be  invoked  to  explain  the  difference  between  the  present  thermal 
attachment  frequency  and  that  measured  by  van  Lint.  The  fractional 
concentration  of  water  vapor  in  the  mixture  is  large  enough  to  ensure 
that  the  energy  exchange  collision  frequency  (which  is  dominated  at  low 
E/N  by  the  water  vapor)  is  sufficiently  larger  than  the  attachment 
frequency  to  prevent  distortion  of  the  electron  energy  distribution  by 
inadequate  thermalization.  No  evidence  for  a  peak  in  the  attachment 
frequency  is  found  in  the  present  work  in  humid  air  in  contrast  to  the 
observations  in  dry  air  (Refs.  28,  29). 

Previous  values  reported  for  the  thermal  attachment  rate 
coefficients  k^  and  k2  are  subject  to  smaller  uncertainties  than  that 
for  k3<  Thus,  the  present  value  of  ka  ■  1.53  x  10-43  m^/s  for  the 
effective  thermal  attachment  rate  coefficient  in  the  humid  air  mixture 
may  be  used  with  the  previous  measurements  of  k^  and  k2  (Ref.  27)  to 
determine  k^  from  Equation  (53).  A  value  of  k-j  ■  (1.0  ±  0.2)  x  10~4^ 
m^/s  is  found  which  is  30%  lower  than  that  reported  previously  (Ref.  2). 
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Further  measurements  in  humid  air  having  higher  water  vapor 
concentrations  would  be  desirable  to  further  reduce  the  uncertainty  in 
the  value  of  kj. 

5.  TWO-BODY  ATTACHMENT  AND  IONIZATION  COEFFICIENTS 

In  paragraph  2  it  has  been  noted  that  the  primary  negative  ions 
formed  in  the  range  of  E/N  >  6  x  10  V *m  are  0  and  H  produced  in 
the  two-body  dissociative  attachment  reactions  described  by  Equations 
(19)  and  (20).  Moreover,  it  appears  that  the  collisional  detachment 
channel  described  by  Equation  (26)  is  the  most  probable  reaction  path  of 
the  H~  ions.  On  the  other  hand,  in  the  case  of  the  0  ions,  collisional 
three-body  conversion  reaction  rates  are  sufficiently  large  that,  over 
an  experimentally  accessible  density  range,  reaction  channels  leading  to 
the  formation  of  the  more  stable  negative  ions  are  favored  over  the 
detachment  reaction  channel. 

For  dry  air  it  has  been-  found  that  the  two-body  attachment  rate 

increases  initially  with  increasing  E/N  over  the  range  60  <  E/N  <  100  x 
—21  2 

10  V'n  but  decreases  rapidly  with  further  increasing  E/N  above  100  x 
-21  2 

10  V*m  (Ref.  1).  In  addition,  it  has  been  noted  that  this  behavior 
is  contrary  to  the  monotonically  increasing  attachment  rate  with 
increasing  E/N  predicted  over  this  range  of  E/N  from  a  Boltzmann 
analysis  using  presently  known  electron  collision  cross  sections  for  N2 
and  O2.  The  increasing  divergence  between  the  measured  and  predicted 
attachment  rates  for  E/N  >  100  x  10”^  V*m^  has  been  attributed  to  the 
increasing  importance  of  collisional  detachment  over  this  range  of  E/N 
(Ref.  1). 

The  presence  of  H2O  in  the  humid  air  mixtures  used  in  the 
present  work  provides  an  additional  important  ion  conversion  channel 
through  the  reaction  described  by  Equation  (25).  Thus,  in  the  case  of 
humid  air,  the  predominant  reaction  channels  for  the  0-  ions  are 


L.  E.  Kline,  private  communication. 
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(21) 


0  +  02  +  M  +  03  +  M 

03"  +  H20  +  M  ♦  O3-  •  H20  +  M  (23) 

0“  +  H20  +  M  +  0“  •  H20  +  M  (25) 

0”+M+0+M+e  (55) 

Further  three-body  reactions  involving  0  »H20  and  O3  *H20  with  H20  lead 
to  higher  order  cluster  ions.  The  formation  of  cluster  ions  decreases 
the  probability  of  collisional  detachment  as  a  result  of  (a)  the 
decreased  mean  energy  of  the  more  massive  cluster  ions  and  (b)  the 
additional  degrees  of  freedom  for  energy  conservation  in  ion-molecule 
collisions  provided  by  successive  detachment  of  H20  molecules. 
Nevertheless,  as  in  the  case  of  dry  air  (Ref.  1),  the  measured 
attachment  coefficient  is  regarded  as  an  effective  coefficient  which  is 
a  function  of  the  true  attachment  coefficient,  the  charge  transfer 
coefficients  and  the  detachment  coefficient.  Increasing  the  gas  density 
favors  the  three-body  ion  conversion  reactions  relative  to  the  two-body 
detachment  reaction.  Thus,  to  minimize  detachment  effects  in  the 
present  study,  measurements  have  been  made  at  the  maximum  values  of  gas 
density  consistent  with  electrical  breakdown  constraints. 

The  values  of  the  net  ionization  coefficient  determined  as  a 

function  of  E/N  from  either  the  negative-ion  or  positive-ion  waveforms 

are  given  in  Table  6  together  with  their  estimated  uncertainties  and 

plotted  as  absolute  values  in  Figure  13  where  they  are  compared  with 

previous  measurements  in  dry  air  (Ref.  1).  For  values  of  E/N  <  70  x 
-21  2 

10  V*m  ,  no  positive  ions  are  detected  and  the  data  in  this  region 
correspond  to  pure  attachment  coefficients. 

The  present  measurements  of  a/n,  which  are  given  in  Table  6, 
have  been  determined  directly  from  the  ratio  of  the  positive  charge  to 
negative  charge  collected  as  described  previously  (Ref.  1),  i.e., 
neglecting  detachment  effects.  The  estimated  uncertainty  in  this  ratio 


TABLE  6.  SUMMARY  OF  THE  SWARM  PARAMETERS  DETERMINED  IN  HUMID  AIR  (2.00%  H2O) 


*/»  , 

10-21V**2 


|4. 44x10 
5.23*10 
3.99*10 
1.80*10 
7.3S*  10 
9.26*  10 
l.U*  10 
1.31*10 
1.40*  10 
|l.7S  *10 
1.79*10 
1.92*10 
2.02  *10 
2.03*10 
2.06  x  10 
2.03*  10 


1.03 

1.03 

1.06 

1.06 

1.03 

1.06 

1.03 

1.03 

• 

1.13 

i.ii 

! 

1.05 

1.06 

1.06 

1.06 

1.16 

1.17 

1.16 

1.16 

1.17 

1.13 

l.tt 

2.47 

2.37 

2.64 

2.21 

2.61 

|  2.32 

1 

2.73 

|  2.31 

2.71 

2.42 

2.79 

2.43 

2.10 

2. a 

2.7* 

2.47 

2.73 

2.40 

2.63 

2.33 

2.23 

2.21 

2.13 

-  0.071  ♦  71 

0.071*  72 

0.37*  12 

-0.142*  3X 

0.142*32 

1.20*  12 

-0.142*  31 

0.31 

0.204.  IX 

2.0  ♦  72 

0. 014.41 

0.61*92 

-0.122*  7X 

0.32 

0.234*102 

2.1  *  UI 

0.132*113 

1.44*12* 

0.40  *32 

2.10 

0.31  *92 

3.3  *112 

0.76  »72 

11. 1  »91 

0.74  ♦  32 

2. *7 

0.31  *62 

3.9  ♦  12 

1.12  i*x 

L7 . 4  ^6X 

1.20  ♦  3X 

4.1 

0.39  *62 

6.4  ♦  12 

i.j*  **x 

24.7  — 1 

2.71  ♦  3X 

7.3 

0.42  *32 

4.2  ♦  71 

j.u 

41  i5I 

4.1  ♦« 

12.1 

0.43  *32 
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—  Dry  Air  ( Davies ) 


Figure  13. 


E/N,  10  “  V-  m 


Comparison  of  the  present  values  (solid  points)  of  the  net 
ionization  coefficient  as  a  function  of  E/N  with  previous 
measurements  (solid  curve)  in  dry  air. 


is  ±  2%.  Combining  these  measurements  with  the  values  of  the  net 
ionization  coefficient  gives  the  coefficients  for  two-body  attachment 
and  ionization  given  in  Table  6  and  Figure  14  where  they  are  compared 

with  previous  measurements  in  dry  air  (Ref.  1).  The  data  for  the 

—  21  2 

individual  coefficients  are  confined  to  values  of  E/N  <  200  x  10  V*m 
for  which  the  ratio  a/n  may  be  determined  accurately.  For  higher  values 
of  E/N,  the  data  correspond  to  net  ionization  coefficients.  Combining 
the  values  of  a/N  and  i/^/N  with  the  values  of  wg  given  in  Table  2 
gives  the  rate  coefficients  for  two-body  attachment  and  ionization  shown 
in  Figure  15  and  in  Tables  7  and  8. 

The  measurements  of  two-body  attachment  in  humid  air  shown  in 

Figures  14  and  15  are  significantly  different  from  those  in  dry  air 

determined  previously  (Ref.  1),  particularly  for  values  of  E/N  > 

-21  2 

100  x  10  V »m  .  However,  the  data  for  humid  air  are  in  good  agreement 
with  the  values  for  0  production  predicted  from  solutions  to  the 

h 

Boltzmann  transport  equation.  This  gives  greater  credance  to  the 
previous  interpretation  that  collisional  detachment  in  dry  air  is 
responsible  for  the  comparatively  lower  values  of  attachment  measured  in 
dry  air.  Although  detachment  effects  are  clearly  much  smaller  in  the 
present  measurements  in  humid  air,  the  measurements  are  still  regarded 
as  effective  values.  For  this  reason,  the  values  of  gas  density  at 
which  the  measurements  were  made  have  been  included  in  Tables  7  and  8. 

Further  measurements  in  humid  air  mixtures  having  different 
water  vapor  concentrations  would  provide  further  information  on  the 
branching  ratio  of  ion  conversion  to  detachment.  It  is  evident  from  the 
previous  measurements  in  dry  air  and  the  present  ones  in  humid  air  that 
this  ratio  is  sensitive  to  variations  in  water  vapor  concentration  over 
the  range  from  0  to  2X.  This  suggests  that  a  series  of  measurements  of 
the  effective  attachment  rate  as  a  function  of  water  vapor  concentration 
at  constant  E/N  could  provide  a  convenient  method  of  determining  the 
branching  ratio  of  ion  conversion  to  detachment  as  a  function  of  E/N. 


L.  E.  Kline,  private  communication. 
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a  (N  / N),  rj  < N  /N),m  *  (N  =  2.688x  10 


E (Ng/N),  104V/m  (N2  =1828x10"" m  ",  226. 5  K,  30 km ) 


E (N, /N),  105 V/ m(N,  =4049xl0”m  ",  216.7  K,15km) 


E  (N„/N),  10  V/m  (N„  =2.688x10  m  ",  273  K,  Okm) 


Humid  Air 

Dry  Air  ( Davies) 


E/N,  10  “  V- 


Figure  15.  Comparison  of  the  present  measurements  of  the  two-body 

attachment  rate  coefficient  (solid  points)  and  ionisation 
rate  coefficient  (open  points)  as  a  function  of  E/N  with 
previous  measurements  (solid  and  dashed  curves,  respectively) 
in  dry  air. 
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v/N.  m  /s 


TABLE  7.  VALUES  OF  THE  TWO-BODY  ATTACHMENT  RATE  IN  HUMID  AIR  £2.00%  h„0) 


N 

1022m"3 

E/N 

10-21V*m2 

m 

v  (2)n  /n 
a  o 

107/s 

Va(2)/N 

10"18m3/s 

138.5 

60.0 

161.3 

1.53 

+ 

8% 

0.57 

+  8% 

138.5 

70.0 

188.2 

3.2 

+ 

6% 

1.20 

+  6% 

138.5 

80.0 

215.0 

5.4 

+ 

7% 

2.0 

+  7% 

76.2 

90.0 

241.9 

7.5 

+ 

11% 

2.8 

+  11% 

32.6 

125.0 

336.0 

14.2 

+ 

11% 

5.3 

+  11% 

32.6 

137.5 

369.6 

15.9 

+ 

8% 

5.9 

+  8% 

32.6 

150.0 

403 

17.2 

+ 

8% 

6.4 

1  + 

00 

7.81 

175.0 

470 

22.0 

+ 

7% 

8.2 

+  7% 

7.81 

200.0 

538 

25.3 

+ 

7% 

9.4 

+  7% 

N  =  2.688  x  10m  J 
o 

Estimated  uncertainties  in  the  measured  quantities  are 
E/N:  +  0.2%;  f :  +  2% 

v  (2)/N,  v  <2>N  /N:  As  noted  in  table 


TABLE  8.  VALUES  OF  THE  IONIZATION  RATE  IN  HUMID  AIR  (2.00%  H20) 


N 

in22  -3 

10  m 

E/N 

.  n— 2L~  2 

10  V’m 

EN0/N 

104V/m 

ViN0/N 

107/s 

Vi/N 

10-18*3/ 

s 

138.5 

80.0 

215.0 

1.64 

+ 

9% 

0.61  + 

9% 

76.2 

90.0 

241.9 

3.87 

+ 

12% 

1.44 

+ 

12% 

32.6 

125.0 

336.0 

29.8 

+ 

9% 

11.1 

+ 

9% 

32.6 

137.5 

369.6 

47 

+ 

6% 

17.5 

+ 

6% 

32.6 

150.0 

403 

72 

+ 

6% 

26.7 

+ 

6% 

7.81 

175.0 

470 

164 

+ 

5% 

61 

+ 

5% 

7.81 

200.0 

538 

306 

+ 

5% 

114 

+ 

5% 

7.81 

225.0 

605 

516 

+ 

5% 

t  192 

+ 

5%. 

3.53 

250.0 

672 

758 

+ 

5% 

t  282 

+ 

5% 

3.53 

275.0 

739 

992 

+ 

5% 

t  369 

+ 

5% 

3.53 

300.0 

806 

1390 

+ 

5% 

t  517 

+ 

5% 

N  -  2.688  x  10~nf3 
o 


J  Net  ionization  rate  coefficient 


Estimated  uncertainties  in  the  measured  quantities  are: 

E/N  ■  +  0.2%;  f  =  +  2% 

v, /N,  v.N  /N:  As  noted  in  table 
i  i  o 


The  analysis  of  such  data  is  particularly  straightforward  (and  less 

ambiguous)  if  the  production  rate  of  0”  is  independent  of  water  vapor 

concentration  at  constant  E/N.  For  this  to  be  true  the  contribution  of 

the  electron  collision  cross  sections  in  water  vapor  in  determining  the 

electron  energy  distribution  in  humid  air  over  the  range  of  E/N  >  60  x 
-21  2 

10  V*m  must  be  small-.  That  such  is  the  case  is  suggested  by  the 

fact  that  the  electron  drift  velocity  in  this  range  of  E/N  is  the  same 

for  dry  air  and  for  humid  air  with  27.  1^0  (cf.  Table  2  and  Figure  9)  and 

(21 

also  by  the  independence  of  values  of  a-rr  ,  determined  by  a  different 
technique  (Ref.  30),  to  water  vapor  concentration  over  the  range  from 
0.06Z  to  1.3%. 

\ 

6.  DISCUSSION 

The  present  data  represent  a  comprehensive  set  of  measurements 
of  the  swarm  parameters  in  humid  air  containing  27.  water  vapor.  This 
concentration  corresponds  to  saturated  air  at  a  temperature  of  291  K. 

The  complete  data  base  for  this  mixture  determined  in  the  present  study 
is  summarized  in  Table  6.  A  plot  of  the  three-  and  two-body  attachment 
frequencies  together  with  the  ionization  frequency  corresponding  to 
atmospheric  density  at  0  km  and  273  K  (2.688  x  10^  m-^)  is  shown  in 
Figure  16  compared  with  previous  measurements  in  dry  air  (Ref.  1). 

Figures  16  and  17  clearly  Identify  the  two  regions  of  E/N  over 

which  the  swarm  parameters  in  air  are  particularly  sensitive  to  water 

vapor  concentration,  i.e.,  E/N  <  25  x  10  ^  y»m^  and  E/N  >  100  x  10  ^ 

2 

V*m  .  It  has  already  been  pointed  out  in  paragraph  5  that  further 
measurements  in  air  at  high  E/N  for  different  water  vapor  concentrations 
would  enable  the  branching  ratio  of  ion  conversion  to  detachment  to  be 
determined  as  a  function  of  E/N.  At  low  values  of  E/N,  additional 
measurements  for  different  water  vapor  concentrations  would  enable  the 
thermal  mobility  of  electrons  and  the  thermal  rate  for  three-body 
attachment  in  dry  air  to  be  determined. 
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iNn  =2.688x  10 


7.  ESTIMATES  OF  MEASUREMENT  UNCERTAINTIES 

The  estimation  of  uncertainties  in  the  present  measurements 
follows  the  detailed  analysis  given  previously  (Ref.  I).  However,  in 
the  case  of  humid  air  there  is  an  additional  uncertainty  introduced  in 
the  measured  coefficients  due  to  the  uncertainty  in  mixture  composition, 
specifically,  the  water  vapor  concentration.  Clearly,  the  effect  of  the 
uncertainty  in  water  vapor  concentration  is  only  significant  under 
conditions  where  the  water  vapor  has  a  significant  effect  on  the 
measured  coefficients.  Thus,  the  estimated- uncertainties  are  the  same 
as  those  for  dry  air  where  there  is  little  or  no  effect  on  the  measured 
coefficients  due  to  the  presence  of  water  vapor.  For  conditions  where 
the  water  vapor  affects  the  coefficients,  estimation  of  the  additional 
uncertainty  requires  a  knowledge  of  the  functional  dependence  of  the 
particular  coefficient  on  water  vapor  concentration  at  each  value  of 
E/N.  In  the  absence  of  such  information,  a  linear  dependence  of  each 
coefficient  on  water  vapor  concentration  has  been  assumed  defined  at 
each  value  of  E/N  by  the  values  determined  in  dry  air  and  in  the  humid 
air  (.21  H20)  mixture.  The  resulting  uncertainties  in  the  swarm 
parameters  for  humid  air  are  estimated  to  be 

UgN ,  UeN/N0:  ±  1 1  for  10"20  <  E/N  <  10-19  V* m2 

±  2%  for  E/N  <  10"2P  V  *m2  and  E/N  >  1CT19  V*m2 
U_N,  U_N/Nq:  ±  1Z 

v.(3)/N2,  v  3N2/N2s  ±  71 

d  d  o 

(a-n)/N, v3^2Vn,  'v^/N:  As  noted  in  Tables  6,  7  and  8 
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IV.  RECOMMENDATIONS 


The  present  measurements  in  humid  air  (2%  1^0)  together  with  the 
previous  data  in  dry  air  (Ref.  1)  provide  comprehensive  sets  of  measure¬ 
ments  for  two  values  of  humidity  of  interest  in  air  chemistry  codes. 
Additional  measurements  in  humid  air  with  different  concentrations  of 
water  vapor  would  be  desirable  to  determine  the  functional  dependence 

of  the  swarm  parameters  on  water  vapor  concentration  over  the  regions 
-21  2  —21  2 

of  E/N  £  25  x  10  V*m  and  >1  100  x  10  V*m  where  the  effects  of  water 
vapor  are  most  significant.  Such  measurements  would  provide  data  on 

(a)  The  mobility  of  thermal  electrons  in  dry  air. 

(b)  The  attachment  rate  of  thermal  electrons  in  dry  air. 

(c)  The  effectiveness  of  water  vapor  as  a  third  body  in  the 
three-body  attachment  process  as  a  function  of  E/N. 

(d)  The  branching  ratio  of  ion  conversion  to  detachment  for 
0  ions  as  a  function  of  E/N. 

All  of  these  processes  may  be  investigated  using  the  present 
drift-tube  facility  which  has  the  capability  for  in-situ  monitoring  of 
the  water  vapor  concentration.  Studies  of  electron  detachment  would  be 
particularly  desirable  in  view  of  the  large  spread  (three  orders  of 
magnitude)  in  the  detachment  rates  that  have  been  reported  in  the 
literature  (Refs.  30-37)  . 

In  addition,  further  measurements  would  be  desirable  on  the 
following  processes: 

(e)  The  effect  on  the  measured  swarm  parameters  of  high-energy 
irradiation  of  the  air  samples  by  cascading  of  the  high- 
energy  electrons  to  energies  in  equilibrium  with  the  drift 
field. 


59 


(f)  The  positive-ion  mobility  at  low  values  of  E/N  below  the  onset 
of  avalanche  ionization. 

(g)  The  recombination  rate  as  a  function  of  E/N. 

These  processes  may  be  investigated  using  the  present  drift-tube 
facility  by  coupling  an  electron  gun  to  the  system  and  injecting  a  high- 
energy  electron  beam  into  the  drift  region.  To  provide  the  flexibility 
in  the  properties  of  the  electron  beam  necessary  for  studying  processes 
(e)-(g)  requires  the  construction  of  an  electron  gun.  The  preliminary 
design  of  a  suitable  gun  has  been  described  previously  (Ref.  38). 
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